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ABSTRACT 



We report the discovery of four dusty cometary tails around low mass stars in 
two young clusters belonging to the W5 star forming region. Fits to the observed 
emission profiles from 24 /im observations with the Spitzer Space Telescope give 
tail lifetimes < 30 Myr, but more likely < 5 Myr. This result suggests that the 
cometary phase is a short lived phenomenon, occurring after photoevaporation 
by a nea rby O star has rernoved gas from the outer disk of a young low mass star 
(see also ISalog et alibood : ISaloe et al.lboosh . 



Subject headings: Circumstellar matter 
winds, outflows — infrared: stars 



stars: pre-main sequence — stars: 



1. Introduction 

Recent observations with the Multiband Imaging Photometer for Spitzer (MIPS) have 
discover ed dusty nebulae a ssociated with several young stars in three Galactic star forming 
regions (IBalog et al.ll2006l ). Appearing as extended, cometary objects pointing away from 
nearby O stars at 24 nm, (and sonietimes at 8 urn), the y resemble the propl yds seen in 
Orion with HST jO'Dell et all Il993l : lO'Dell fc W^ll994j ). ISalog et all tood ) show that 
these cometary tails are probably dust swept out of a young circumstellar disk by an O star. 



Balog et al.l (120061 ) show that the emission profiles at 24 and 8 /im can be explained by 



a model where gas is photoevaporated from a young protoplanetary disk. As the gas leaves, 
entrained dust also escapes and is then blown away from the disk by radiation pressure from 
nearby O stars. The dust is swept into a cometary morphology which glows in the mid-IR 
thermally as it reprocesses t he inc i dent UV radiation to longer wavelengths. iBalog et a 



and 



fl2006h cite IJohnstone et al.l (119981) . iRichhng fc Yorkd (l2000h . iHoUenbach fc Adamd (120041 ) 



Matsuyama et al.l (120031 ) . who find that the portion of a protoplanetary disk likely 
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responsible for emission at 24 /im 
yr. 



^2-50 AU from the star) is removed on timescales ~3 x 10^ 



We have recently discovered four d usty cometar y struc tures in the W5 star forming 
region. Here we expand on the work of iBalog et al.l (120061 ) and place constraints on the 
lifetimes and physical nature of these interesting objects. 



Observations 



Koenig et al.l (120081 ) describe our observations of W5. We observed a 2°xl.5° field with 
the Infrared Array Camera (IRAC) and MIPS in 2006 and 2007. Figure 1 shows the MIPS 
24 /im image mosaic. Images were taken in scan map mode using the medium scan speed 
for an average exposure time of 41.9 seconds per pixel for the combined frames. The raw 
MIPS data were processed with pipeline version S13.2.0. The final mosaics were produced 
using the MIPS instrument team Data Analysis Tool, which calibrates the data and applie s 
a distortion correction to each individual exposure bef ore combining (I Gordon et al.ll2005l ). 
We used the dustergrinder software tools (described in iGutermuth et al.l l2008l ) to produce 
final image mosaics from the IRAC data in each of the four wavelength bands. 

Figure 2 shows the 24 /im cometary objects at larger scale. At wavelengths shorter than 
24 iim, these objects are point-like or not detected. At 70 and 160 /im, these objects are 
also invisible. Table 1 summarises the properties of the four objects in W5. The quoted 24 
lim flux is for the heads of the cometary features alone. 



3. Analysis 



We follow iBalog et al. I J2006h and analyze these objects using an adaptation of the Vega 



debris disk model of ISu et al.l (j2005l ). We assume the outward flow has roughly constant 
radial density: n{r) ~ constant. We fit the emission profile of a 5 pixel (6.23") wide strip 
along the symmetric axis of the tail. To estimate and subtract the background, we extend 
the strip to reach regions free of diffuse emission and point sources. 

The fit to the resultant surface brightness profile is based on optically thin emission 
from dust at its equilibrium temperature. The emission depends on dust composition, dust 
grain size, and the angle that the tail makes with the plane of the sky (inclination), where 
0° corresponds to a tail in the plane of the sky. We adopt astronomical silicates for the dust 



composition. Due to the degeneracy between tail inchnation and dust grain size (IBalog et al. 



20061 ). and the lack of reliable constraints on these parameters, we assume an inclination of 
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0° and fit the profile adopting a single grain size. 

The fits to the 24 pm emission profiles are shown in Figure 3. Error bars on the profile 
points are derived by extracting the emission profile to either side of the main axis and 
taking the difference between the two extractions as the error at each point. Tails A, B and 
C are formed by the multiple O star system HD 17505 which contains at least four O stars 



( iHillwig et al.l 120061 ). We model the radiation output of this system as a single star with 
effective temperature T^ff = 35000 K and radius R = 21Rq. Tail D is most likely formed 
by the 07V star HD 18326. We adopt T^ff = 36000 K and i? = 9ARq. The model provides 
reasonable fits to the dust emission profiles at 24 /um, however all of the fits diverge from 
the data at either end due to contamination from nearby stars. 

At 8 iJ,m, our model predicts a flux of 0.0025 mJy for tail A. This flux is consistent with 
the upper limit listed in Table 1. For tails B and D the model underpredicts the 8 /xm flux by 
more than an order of magnitude. Tail C is a more complex case. Due to crowding, we could 
not derive a precise flux for this object. The model predicts a flux >5 times the upper limit 
quoted in Table 1 for tail C. For tails B and D, the enhanced level of 8 /im emission in our 
da ta is probably prod uced by stochastically heated small grains, ~ loA in size as suggested 



by iBalog et al.l (120061 ). To test whether this is also true for tails A and C requires improved 



measurements of the flux at 8 fim. 

Optical and near infrared observations of the dusty tails in NGC 2264, IC 1396 and 
NGC 2244 (Balog et al. 2008) suggest that these tail structures are strongly deflcient in gas. 
If gas does not entrain the dust away from the circumstellar disk, another mechanism must 
produce these structures. We propose that radiation pressure from nearby O stars ejects the 
dust directly. To test this idea, we calculate the radius, Vevac within the stellar disk where 
radiatio n pressure exceeds the gravitational binding energy or dust in orbit around its host 
star (see 



Burns et al.lll979l ). For a star of mass nistar- 



'SGpnd? c a nistar ^^-^ 



3LoQ 



pr 



where d is the distance from the tail head to an O star with luminosity Lq, and Qpr is the 
frequency averaged radiation pressure coefficient {Qpr ~1 for an O star and astronomical 
silicate dust grains; Laor & Draine 1993)0. For a 1 M0 host star and grains of size a = 0.01 
fim, and density p = 3.3 g cm~'^, Vgyac = 9-35 AU. For 10 pm grains, Vevac ^ 500 AU, the 



^The expression for r^vac assumes grains produced from collisions of larger objects. For primordial grains 
orbiting the central star, r^vac is a factor of y/2 larger. 
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typical protostellar disk radius (lAndrews fc Williamsl 120071 ) . Thus radiation pressure can 
eject grains with a <3 /xm. 

Table 2 hsts two values for the tail mass. The lower limit assumes a single grain size of 
0.01 fim. The upper limit assumes a size distribution fr om 0.01 /i r n up t o th e grain size at 
which 

fevac ~ 350 AU. We adopt the size distribution of iDohnany j (119691 ) and I .anaka et al.l 



(1l996l ). N{a) oc a^'^'^ for a grain population formed by a coUisional fragmentation cascade 
in a protoplanetary disk, normalize the mass in 0.1 /xm grains require d to fit the s urface 
brig htness distribution (F i g. 3) and integrate over grain size. Recently, IJorda et al.l (120071 ) 
and iHolsapple fc HousenI (120071 ) found a different power law (—3.1) for dust grains in the 
Deep Impact ejecta experiment, which would result in a factor ~10 greater tail mass. 

To derive a simple mass-loss rate for each tail, we divide the estimated tail mass by 
the projected tail length over the flow velocity (follo wing Balog et al. 2008 we use the 
drift velocity given by eqn. B.4 of lAbergel et al.l 120021). To derive an ex pect ed lifetime, we 



adopt a minimum mass solar nebula as specified by I Weidenschillingj (119771 ) and lHayashi et al. 
( 1l985l ) as a dust reservoir and find a total mass Mtot = 0.034 Mq between radii 20-350 AU. 



Assuming the dust mass is a factor 100 smaller, the mas s loss rates yield lifetimes of 0.3-30 
Myr (Table 2). We note that lAndrews &: Williamsl (120071 ) find a range in disk masses around 
young stars, such that Mtot could be as much as a factor 10 larger. 



Discussion 



The initial phase of a young stellar disk's evolution is the optically-thick phase. This 
period lasts a few to ~10 Myr; during this time, excess emission from dust in the circumstellar 
disk in t he Spitzer IRAC bands (3-8 fim ) is visible, as are signatures of accretion onto 
the star (ISicilia-Aguilar et al.l l2005l . |2006| ). The cometary tails in W5 are unlikely to be 
older than 5 Myr. Since massive O stars are present and non-thermal radiation is largely 
absent (which would be the signature of a supernova haying previously occurred ), 5 Myr 
is a conservative age upper limit for W5 (iKarr fc Martini l2003l : IVallee et al.lll979l ). During 
this early phase, UV radiation from nearby O stars will photoevaporat e the outer disk of 



stars in their cluste r s, creating p r oplyd s such as those seen in Orion (lO'Dell et al.l Il993 



O'Dell & Wen 


1994 


)• 


Balog et al. 


(2008) 



micron-size dust grains are replenished by collisions between already-formed planetesimals 
and then blown out by radiation pressure, creating the cometa ry tails seen at 24 irni. We 



consider this scenario somewhat unlikely in such young disks. iKenyon &: Bromley! (1200 ll ) 
find that a size distribution of bodies can stir up planetesimals enough to begin a collisional 
cascade capable of making small dust grains in the outer 30-150 AU of a disk, only if it 
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contains objects of 500 km radius or larger, whi ch need at least ~10 Myr to form at this 
distance from the star (IKenyon fc Bromleyll2008l ). Since we find that the dust grains that 
create the tails come from beyond ~10 AU they must have remained in the disk and survived 
being completely dragged out as the gas departs during photoevaporation. 

We note that some of the objects could be older than the O stars and have migrated 
in to their present location. This scenario requires an earlier generation of stars (5-10 Myr 
ago) to have formed in a cluster of ~1000 stars within ^5 pc of HD 17505 or HD 18326. 
In addition, faster planet formation (e.g. by d isk instability, or induced by the radiative 
environment Boss 2006 : Throop fc Ballyl bood ) might allow the planetesimal stirring and 



coUisional cascade mechanism to operate and create the small dust grains that make up the 
comet ary tails. The dust may not r emain long enough in such a disk to make detection likely 
however (IKenyon fc Bromleyll2008l ). 



Tails A, B and C belong to the large double cluster of young stellar objects (YSO) 
surrounding O stars HD 17505 and HD 17520. T his cluster ha s a pr ojected radius on the 
sky of 3 pc, within which are 146 stars classified by lKoenig et al.l (120081 ) as class II stars with 
disks, i.e. stars in the optically thick disk phase, including tails A, B and C. If the cluster is 
spherical and uniformly distributed we would expect 1.54 objects within a spherical radius 
0.67 pc (the distance of tail A from HD 17505). Similarly tail D belongs to the cluster around 
HD 18326 which has a radius of 1.8 pc and 87 YSO/optically thick disk objects. Within 
a 0.22 pc radius (that would encompass tail D) we would predict 0.16 objects. Given that 
projection effects mean we should see only one third of these objects as extended, it is likely 
that all objects with disks capable of producing cometary tails within 0.7 pc of an O star 
are doing so. 

Not all O stars appear to create tails however. Within W5, 2 O star systems (several 
are binary or multiple) out of 5 possess cometary objects at 24 /im. Spitzer observations of 
star forming r egions with MIPS (GTO program ID 58) have found 3 more out of a sample 
of 20 O stars (iBalog et al.ll2006l ). This sample includes the Orion Trapezium re gion which is 



satur ated at 24 /im but does contain the well known optically visible proplyds (lO'Dell et al. 



19931 ). The timescale of the mid- infrared cometary phase may thus be closer to the short 
end of our estimated range. 



5. Conclusions 



We have discovered four dusty cometary tails in W5 at 24 /im with Spitzer MIPS. 
An apparent deficit of gas suggests an origin in radiation pressure blowout of dust from a 
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young stellar disk by nearby O stars. Future observations with high sensitivity, high spatial 
resolution mid-infrared and sub-millimeter instruments will help constrain their nature. 

Facilities: 2MASS (JHRg), Spitzer (IRAC, MIPS) 
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Table 1. 





R.A. 


Declination 


./ 


H 


Ks 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] 




.T2000.0 


J2000.0 


(m,7y) 


{m,Iy) 


(mJy) 


(mJy) 


{m,Iy) 


(mJy) 


(m,7y) 


(mJy) 


A 


2 51 15.17 


60 24 20.31 








<0.006 


<0.007 


<0.035 


<0.048 


19±3 


B 


2 51 13.53 


60 25 02.36 








0.40±0.04 


0.27±0.02 


0.15±0.05 


0.19±0.05 


16±2 


C 


2 51 10.49 


60 25 07.15 








<0.006 


<0.007 


<0.035 


<0.048 


26±2 


D 


2 59 20.91 


60 34 14.82 


0.43±0.09 


0.54±0.12 


0.69±0.12 


0.42±0.04 


0.32±0.04 


0.18±0.08 


0.59±0.06 


37±3 



Note. — JHKq photometry obtained from 2MASS catalog jSkrutskie et al.|200^ . We obtained the remaining photometry in 6'pi^2:er bands 
at tail head position, using aperture photometry. 



Table 2. 





Tail Length 


Tevac (0.01//m grains) 


o) 

"'max 


Tail Dust Mass"'- 


Vdrift 


Lifetime 




(pc) 


(AU) 




(10-6 Mo) 


km s"^ 


(Myr) 


A 


0.3 


34.5 


1 


0.15 - 1.5 


40 


1.6 - 15.8 


B 


0.49 


20.2 


3 


0.08 - 1.3 


70 


1.6 - 30.1 


C 


0.62 


9.3 


10 


0.06 - 4.7 


150 


0.3 - 24.5 


D 


0.29 


23.6 


3 


0.14 - 2.4 


55 


0.7- 12.6 



Note. — ''^Maximum grain size evacuated from disk under assumption that no significant 
material exists beyond 350 AU. ^Tail dust masses for assumed single grain size and for power- 
law grain size distribution. 
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3:02:24.0 99:12.0 2:56:00.0 52:48.0 49:30.0 46:24.0 

Right ascension 



Fig. 1. — 24 //m mosaic of W5, as imaged by Spitzer/MIPS. Two bubbles are seen in 
projection on the sky. Cometary sources are seen in the vicinity of the primary ionizing O 
stars. Expanded views of the 2 boxed regions are shown in Figure 2. The left hand box 
covers HD 18326, and the right hand box HD 17505. 
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Fig. 2. — Expanded views of the boxed regions in Figure 1 in exaggerated color to show 
the cometary sources (marked with crosses) and ionizing stars (marked by tics). Several 
cometary sources point to the multiple O star system HD 17505 (upper panel, see Hillwig 
et al. 2006). One cometary object is evident near HD 18326 (lower panel). 



- 11 - 




Fig. 3. — Fits to the 24 fim brightness distributions of the four cometary tails in W5. 
Contamination from neighboring sources and surrounding diffuse emission hkely causes much 
of the discrepancy between the data and the fits to tails B, C and D. Data points are shown 
by diamond points with error bars. The best fit is the solid blue line passing through the 
points. 



